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on behalf of the Fermi LAT and GBM collaborations 

The bright short-hard GRB 090510 was observed by both Swift and Fermi telescopes. The study 
of the prompt emission by Fermi revealed an additional high-energy spectral component, the largest 
lower limit ever on the bulk Lorentz factor in a short GRB jet, and brought the most stringent 
constraint ever on linear Lorentz invariance violation models. The fast repoint and follow-up by 
both telescopes allowed the first multiwavelength study of a GRB afterglow from optical range to 
several GeV. This long-lived emission has been studied in the framework of the internal shock and 
external shock models. 



I. OBSERVATIONS 

On May 10th, 2009 at 00:23:00.48 UT, the Burst 
Alert Telescope (BAT) onboard the Swift observatory 
detected a bright Gamma- Ray Burst (GRB) 1], and 
autonomously repointcd towards it within 100 s, al- 
lowing for over 2 days of follow-up observation by the 
Ultra- Violet and Optical Telescope (UVOT) and the 
X-Ray Telescope (XRT). An afterglow was detected 
in ultra-violet and X-ray, which analysis provided a 
location RA = 333.55227, Dec = -26.5827, with an 
error of 1.4 arcsec (90% statistical) @. VLT after- 
glow observations_provided a spectroscopic redshift z 
= 0.903 ± 0.003 Q. 

On the same day at 00:22:59.97 UT, the Gamma- 
ray Burst Monitor (GBM) onboard Fermi detected a 
bright GRB |4j. A small and short spike at trigger 
time was followed 0.5 s later by 0.5 s-long emission 
made of 7 pulses. This GRB also triggered the Fermi 
Large Area Telescope (LAT) onboard algorithm Q 
which located it at RA = 333.4°, Dec = -26.767°, 
with an error of 7 arcmin (statistical), consistent with 
the GBM and Swift detection. The spacecraft au- 
tonomously repointed towards the source for 5 hours 
of follow-up observation. The prompt emission was 
coincident with the GBM main emission. More than 
150 photons above 100 MeV were observed, and more 
than 20 events above 1 GeV in the first minute after 
trigger (see fig. [!} [6| . The follow-up observation al- 
lowed the detection of a significant GeV emission up 
to 200 s after trigger. 

Also AGILE Q, Integral-SPI and Suzaku-WAM 
detected GRB 090510 prompt emission. According to 
all observations GRB 090510 was a very bright short 
GRB (table HI. 

In the following we focus on the prompt emission 
seen by Fermi (section [TTJ) and the multiwavelength 
afterglow study based on Swift and Fermi data (sec- 
tion IIIip . Convention : spectral and temporal in- 
dices are defined according to to the fluence evolution 



TABLE I: Duration (T 90 and T 50 ) of GRB 090510 prompt 
emission observed by different instruments. [l(| 



Instrument 


T go (s) 


T 50 (s) 


GBM/Nal 3,6,7 


2.1 


0.2 


Swift/B AT 


4.0 


0.7 


Integral-SPI 


2.5 


0.1 


Suzaku-WAM 


5.8 


0.5 



description : F oc t 



II. THE 7- RAY PROMPT EMISSION 



A. Spectroscopy 



The spectral analysis of the prompt emission was 
performed with RMfit, combining GBM and LAT 
data (fig.©. 

The time-integrated spectrum shows a significant 
additional high-energy power-law component (N a = 
5.6, (5 PL = 0.62 ± 0.03), which carries 37% of the 
total 7-ray fluence. 

The time-resolved spectroscopy shows the late ar- 
rival of the LAT emission and of the additional com- 
ponent. In bin 'a' (0.5s - 0.6s after Fermi trigger), 
a Band function with steep fiBand > 4 is fitted. In 
bin 'b' (0.6s - 0.8s), a significant additional power-law 
component is found (fiph — 0-66 ± 0.04). In bin 'c' 
(0.8s - 0.9s) this high-energy component can be fitted 
but not significantly. In bin 'd' (0.9s - 1.0s) the LAT 
data only are fitted by a power-law {fipL = 0.9 ±0.2). 
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TABLE II: To , 



values from GRB 090510 observations. 



0.5 1 1.5 

Time since GBM trigger (May 10, 2009, 00:22:59.97 UT) (s) 



FIG. 1: GRB 090510 observation by Fermi [<|. 
Panel (a) : LAT events passing the onboard (blue) or on- 
ground (red) photon selections. The lines show the arrival 
time dependence on photon energy for the events asso- 
ciated to the 31 GeV photon, according to linear (solid 
lines) or quadratic (dashed lines) energy dependence in 
the LIV model (see section |IID[) . Different parts of the 
lower-energy emission are considered : full GBM emis- 
sion (black) , bulk GBM emission (red) , bulk LAT emission 
(green), coicident GBM peak (blue). 

Panels (b — f) : multi-instrument lightcurve. Curves 
(b - d) are background-subtracted. Curve (d) shows all 
LAT events passing the onboard photon selection. Curves 
(e - f) show on-ground selected photons inside the region 
of interest, used for the spectral analysis. Panel (f) also 
shows event energies. 



B. Limits on the jet's initial bulk Lorentz factor 

To 



The observation of high-energy events allows to con- 
strain the jet's initial bulk Lorentz factor Tq. Low- 
and high-energy photons emitted in the same physi- 
cal region of the jet can interact to produce pairs. 
The opacity to pair-production is strongly reduced if 
the emission region moves towards the observer with a 
high bulk Lorentz factor. The optical depth to e^-pair 
production at a given photon energy r 77 (i?) depends 



Time bin 


Spectral model 


t v (ms) 


E max (GeV) 




b 


Band + PL 


14 ± 2 


3.43 


951 ± 38 


d 


Band 


12 ± 2 


30.5 


1324 ± 50 


d 


Band + PL 


12 ± 2 


30.5 


1218 ± 61 



on the 7-ray emission spectral shape and variability 
time scale t v , the GRB redshift and the jet's bulk 
Lorentz factor. The condition r 77 (i? < E max ) < 1, 
where E max is the maximal observed photon energy, 
yields a lower limit on Tq. 

For different epochs of the LAT emission, estimates 
of t v and spectral models lead to several possible lower 
limits, all of order T min ~ 10 3 (see table UTJ) . This 
is the highest lower limit ever set on a GRB, and the 
most rapid ejection ever observed from a short GRB. 



C. Interpretation of the high-energy spectrum 

The physical interpretation of the high-energy 
prompt spectrum of GRB 090510 is described in de- 
tail in [lOj . The lower-energy part of the prompt emis- 
sion spectrum (< 10 MeV) can be explained by non- 
thermal synchrotron processes, as well as photospheric 
thermal emission. 

In a non-thermal leptonic model framework, Syn- 
chrotron Self-Compton (SSC) emission is expected at 
GeV-TeV energies and can result in the observed ad- 
ditional power-law component. Detailed simulations 
of SSC emission have been performed, using the ob- 
served Band component as the seed population, taking 
into account internal 77 opacity, self- absorption, ra- 
diative escapes, and considering either Thomson or 
Klein-Nishina scattering regimes. A low magnetic 
field (B « 10 3 G) favors Compton scattering and 
can explain the observed bright and hard high-energy 
component. It would also slow down the cascade 
formation which could explain the late onset of the 
high-energy emission. The interpretation of the high- 
energy emission as forward shock emission from an 
early afterglow was also considered to explain the de- 
layed onset of the high-energy emission. However, 
this scenario assumes a very rapid ejection (T ~ 2000 
- 4000), and implies a high-energy component spec- 
tral index of (3pl — 1. This spectrum is consistent 
with the long-lived spectrum, but not with the addi- 
tional component found in bin 'b' (/3pl = 0.66). In 
both cases, the hard spectral index of the soft photons 
(oiBand — —0.48) is not explained by the standard syn- 
chrotron mechanism. 

In hadronic models, secondary e pairs come from 
the decay of pions formed through photohadronic pro- 
cesses and by the 7—7 attenuation of synchrotron pho- 
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Time- integrated photon spectrum (0.5 s - 1 .0 s) 
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FIG. 2: GRB 090510 prompt spectrum 0. 
Panel (b) : time-integrated Spectral Energy Distribution 
(SED). The Band function and power-law components are 
shown. Their sum is drawn in thick line with 68% con- 
fidence level contours. Panel (c) : time-resolved SED 
with 68% confidence level contours. In bin 'a' no signal is 
detected in the LAT. Bins 'b' and 'c' show an additional 
power-law component. No signal is detected in the GBM 
after 0.9 s (bin 'd'). 



tons radiated by protons and ions. Synchrotron radi- 
ation or inverse Compton scattering from these pairs 
can produce the observed high-energy 7-ray emission. 
In the case of photohadronic models, a stronger mag- 
netic field makes secondary pion production more ef- 
ficient through more rapid hadronic acceleration, but 
results in a high-energy component spectral index 
(3 pl ~ 1, significantly softer than observed. On the 
other hand, inverse Compton scattering of secondary 
pairs can produce a high-energy component as hard as 
the one observed, but this scenario requires a very high 
isotropic-equivalent power Pi SO ~ 10 55 er<7.s _1 . Pro- 
ton synchrotron models also require a high isotropic 
equivalent energy, and a high magnetic field. How- 
ever, the apparent Ei SO problem can be solved if the 
jet is strongly collimated or if the true Lorentz fac- 
tor r is a factor ~2 less than the minimum values 
ro.rmn estimated through the 7 — 7 opacity constraint. 
If r > ro jm i n , then a proton synchrotron model is 
strongly disfavored, though note that GRBs with a 
very high isotropic energy have been observed. 



D. Limits on Lorentz Invariance Violation (LIV) 

Some quantum gravity models allow for a depen- 
dence of the photons speed v p h on their energy, i.e. 
their arrival time can be written as a development in 
their energy. In this framework, a low-energy (Eg) and 
a high-energy (Eh) photons emitted together would 



TABLE III: Limits on linear LIV derived from GRB 
090510 observations 0. 



\At/AE\ or |At| 


Mqg,i/M P i 


inch 




Method 


< 30 ms/GeV 


> 1.22 




± 1 


LAT lag analysis 


< 859 ms 


> 1.19 




1 


not emitted 










before GBM onset 


< 10 ms 


> 102 




± 1 


GBM pulse width 



arrive at different times, with the dominant LIV de- 
lay term : 



At = 



(1+n) (El-E^p 
2H Q (M QG ^Y 
(l + z') n 

^n m (i + z'f Wl 



dz' 



where z is the GRB redshift, and n = 1 or 2. A stan- 
dard cosmology (h, fi ro , fi A ) = (0.71, 0.27, 0.73) has 
been assumed for the computations. The difference 
in arrival times can be of any sign s n , depending on 
the models. The observation of high-energy photons 
from distant sources allows to put constraints on the 
energy scales Mqq n that appear in the development. 

GRB 090510 has a redshift of z = 0.903 ± 0.03, 
and its emission in the LAT energy range includes a 
photon of energy E h = 30.5±|;| GeV detected 0.829 s 
after the GBM trigger. This event's topology makes 
it a very good photon candidate. The expected back- 
ground rate and this event's direction proximity to 
GRB 090510 location allowed us to associate this 
event to the burst with high confidence (A^ =4.4 to 
5.6 depending on the selections used). 

The observed lightcurve allows different measure- 
ments of the temporal delay for this event, depending 
on which part of the lower energy emission it is asso- 
ciated to. As described in detail in [9j, several lower 
limits on the linear term's energy scale Mqg,i could be 
derived (see table IIIip . All require Mqq : i > Mpi anc k, 
which is a strong constraint on linear LIV models. 

Lower limits on Mqq^ were also derived, but they 
are much less constraining. For more details on LIV 
limits derivation, see 12 1. 



III. 



AFTERGLOW MULTI- WAVELENGTH 
STUDY 

A. eV to GeV observations 



The fast repoint by Swift and the ARR performed 
by Fermi provided a long-lasting observation of GRB 
090510 afterglow, from the optical range to several 
GeV, interrupted by Earth occultation of the GRB to 
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FIG. 3: GRB 090510 Swift and Fermi observations 
Top panel : LAT flux above 100 MeV and best fit to 
the flux decay (line). Bottom panel : energy flux den- 
sities averaged in the observed energy bands: BAT (15 
keV - 350 keV, stars); XRT (0.2 keV - 10 keV, crosses); 
UVOT renormalised to white (diamonds); LAT (100 MeV 
- 4 GeV, filled squares; the average spectral index was used 
to convert from photon to energy flux) with upper limits 
for P — 1.1 (triangles). The prompt emission is shown 
for comparison: GBM (8 keV - 1 MeV, circles), LAT (100 
MeV 4 GeV, empty squares). XRT lightcurve is obtained 
as in fLU . [Ta |. All other data are shown with 68% error 
bars or 95% confidence level upper limits. 



both observatories. The fluxes observed by all instru- 
ments are reported in fig. El 

UVOT and XRT observations of the afterglow 
started ~ 100 s after the prompt emission trigger, and 
lasted over 2 days. The optical and UV lightcurve is 
best fit by a smoothly broken power-law with an ini- 
tial rise : ao P t,i — — 0.50l n 'iq- The smooth break 



J -o.i3- 

1.58^0 37 k s is followed by a shallow decay 

j+O.'ll 



at topt — 

&Opt,2 = 1-13^qiq. The X-ray lightcurve is best fit 



by a broken power-law, with an initial shallow decay 
a x ,i = 0.74 ± 0.03, a break at t x = 1.43±g;?§ ks and 
a steep late decay ax, 2 = 2.18 ± 0.10. 

The beginning of the temporally-extended emission 
in the LAT has been chosen as the end of the prompt 
emission seen by the GBM (0.9 s after GBM trigger, 
i.e. 0.38 s after BAT trigger). A significant emis- 
sion was detected up to 200s after trigger and up to 
4 GeV. This observation was divided in time bins, in 
all of them an unbinned likelihood spectral analysis 
was performed using diffuse class events The 
flux decay is best fit by a simple power-law of index 



3 

"m O 



E S 



.-t-t 




10- 
10 

1 10- 7 



10- 

10" 1 



10 3 0.01 0.1 1 10 100 1000 10 4 10 5 
Energy (keV) 


10 6 


10 ? 


r 


1 ' """1 1 1 1 '1 ' ' 




: 1 H-fa 












XRT '•»., 






- UVOT 














LAT 




















10" 3 10" 2 10" 1 1 10 10 z 10 3 10" 10 5 10 6 
Energy (keV) 



FIG. 4: GRB 090510 UVOT-XRT-LAT SEDs [TJ]. 
Top : SEDs at different epochs, with the best fit shown 
(P2 and 03 constrained, see text). The butterfly at 100s 
indicates the 68% confidence level region for the LAT flux, 
obtained from an unbinned likelihood analysis (95% error 
bar at 100 MeV is shown) . Successive SEDs in time order 
are rescaled by 1:1, 1:10, 1:100, 1:1000, 1:10000. Bottom 
: SED at 100s, with the best fit shown. The /3 2 - j3 3 
constraint is relaxed. 68% error bars are shown for UVOT 
and XRT data. A 68% confidence level region is drawn in 
the LAT energy range. 



a 7 — 1.38 ± 0.07, with no significant feature or break. 
The power-law spectrum shows no significant evolu- 
tion over time and has an average index j3~ f = 1.1 ±0.1. 



B. Discussion on the afterglow origin 

The interpretation of the long-lasting multi- 
waveband emission fron GRB 090510 is discussed in 
detail in [111 ]. Two scenarios are considered in the 
fireball model frame. 

In the first scenario, X-ray and 7-ray emissions are 
due to internal shocks, and the optical and UV emis- 
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sion to the forward shock. The fluxes measured at 
100 s after trigger in the X-ray and 7-ray ranges are 
consistent with an internal shock origin, with some 
fine tuning required. For instance this scenario as- 
sumes that the initial rise of the optical emission is due 
to the forward shock onset, and therefore is steeper 
than what was observed. But the observation may 
have caught the end of the steep rise phase. This in- 
terpretation also requires a very low ambient density 
n ~ 10 -6 cto -3 , which is low, even for short GRBs. 
The assumed initial bulk Lorentz factor Tq ~ 10 3 is 
consistent with the prompt emission analysis (see sec- 
tion htbji. 

Another possibility is that the full long-lasting emis- 
sion comes from the forward shock region. This 
model predicts a broad spectrum with a doubly- 
broken power-law shape and constraints on the indices 
: fix = —1/3, /?3 = P2 + 1/2. Five successive Spec- 
tral Energy Distributions (SED) have been built for 5 
different observation dates (100 s, 150 s, 700 s, 1000 
s, 1400 s after trigger), including X-ray and UV data, 
as well as LAT data for the SED at 150 s (see fig. 0}. 
A fit of these data yields a good agreement with the 
aforementioned model, with 02 = 0.78 ± 0.04. The 
low-energy break decreases with time, from 0.43 keV 
to <0.01 keV. The high-energy break is poorly con- 
strained (e [10 — 130] McV) but is confirmed by a fit of 
the 100 s SED alone (N a > 4.8) where the constraint 
on 02 and P3 is relaxed but still holds within error 
bars. As a result, the afterglow emission spectrum 
is well described by the forward shock model over 9 
energy decades. The lightcurves observed before 1 ks 
after trigger are also consistent with this model. How- 
ever, the early onset of the forward shock emission 
requires a very high bulk Lorentz factor To > 5800, 
and some temporal emission properties are not well 
explained, e.g. ax,i is too shallow and ao P t,2 7^ ctx,2- 
Theoretical extensions of the model can alleviate these 
problems. E.g., a phase of energy injection or an evo- 
lution of the microphysical parameters of the blast 



wave may cause an early shallow decay of the X-ray 
flux ; also the difference between X-ray and optical 
late decay slopes could be explained by dynamical ef- 
fects. 



IV. CONCLUSION 

GRB 090510 observation was remarkable for sev- 
eral reasons. First, it was the first short GRB of 
known redshift (z = 0.903) observed at GeV energies. 
It yielded the highest energy photon ever observed 
from a short GRB (30.5 GeV). The observation of 
high-energy photons from this GRB allowed to put 
strong and robust constraints on linear LIV models 
(Mqg,i > Mpi anc k required) and on the jet's velocity 

(ro,min ~ 10 3 ). 

This GRB also yielded the first clear evidence of 
an additional high-energy spectral component (N a = 
5.6 on the time- integrated prompt spectrum ), w hich 
was to be observed in several other bursts fill [i~7| . 
This high-energy spectral component is a hint for SSC 
radiation or a possible UHECRs production in GRBs. 

Finally, this energetic short GRB showed a bright 
optical and X-ray afterglow, as well as a long-lived 
GeV emission, which could be observed thanks to 
Swift and Fermi fast repoint abilities. The temporal 
and spectral properties of this long-lasting emission 
could be studied over 9 energy decades and can be 
explained in at least two ways : a combination of in- 
ternal shock and forward shock emission reproduces 
well the fluxes observed although some fine tuning is 
required, and a forward-shock only origin can repro- 
duce the observed spectra although it requires some 
theoretical extensions. Such joint Swift-Fermi obser- 
vations are very promising for understanding GRB af- 
terglow physics. 
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